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Berry Curvature, Triangle Anomalies, and the Chiral Magnetic Effect in Fermi Liquids 
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In a three-dimensional Fermi liquid, quasiparticles near the Fermi surface may possess a Berry 
curvature. We show that if the Berry curvature has a nonvanishing flux through the Fermi sur- 
face, the particle number associated with this Fermi surface has a triangle anomaly in external 
electromagnetic fields. We show how Landau's Fermi liquid theory should be modified to take into 
account the Berry curvature. We show that the "chiral magnetic effect" also emerges from the Berry 
curvature flux. 
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Introduction. — Recently there has been a lot of interest 
in the effect of the Berry phase and Berry curvature on 
the physics of the electron Fermi liquid. The standard 
theory of Fermi liquids, developed by Landau [1], as- 
sumes that the low-energy degrees of freedom in a Fermi 
liquid are the fermion quasiparticles, whose distribution 
function in phase space satisfies a kinetic equation. In 
many cases the semiclassical motion of a wave packet of 
electrons in a crystal should include an extra term due 
to the Berry phase, expressible in terms of the electronic 
Bloch wave functions [2]. Such a term should alter the 
standard kinetics of Fermi liquids (see also below); in- 
cluding this term leads to an interpretation of the anoma- 
lous Hall conductivity in terms of Fermi surface proper- 
ties [3]. 

In this Letter we show the connection between the 
Berry curvature on the Fermi surface and triangle anoma- 
lies. Let us first notice that the total flux of Berry cur- 
vature through a given Fermi surface does not need to 
vanish, but can be a multiple of the flux quantum. One 
possible case is doped Weyl semimetals [4-6], but the 
discussion does not depend on the origin of the Berry 
curvature flux. We will show that if there are k quanta 
of Berry curvature flux through a given Fermi surface, 
then the number of particles associated with this Fermi 
surface (which is proportional to its volume) is not con- 
served in the presence of the external electromagnetic 
field, 



- + V.j^^E.B. 



(1) 



Charge conservation is ensured by the vanishing of the 
sum of fc's of all Fermi surfaces. 

Equation (1) is exactly the equation of triangle anoma- 
lies in relativistic quantum field theory [7, 8]. It is 
therefore expected to hold for a Fermi gas of relativis- 
tic fermions at finite density. Indeed, the Berry curva- 
ture of a relativistic fermion has the form of the field of 
a magnetic monopole in momentum space, and k ^ ±1 
for right- (left-)handed fermions. The statement (1) goes 
further by tying anomalies with Fermi surface properties 
only. In this way, we demonstrate that axial anomalies 
are properties of Fermi liquids with Berry curvature flux. 



even when the original particles interact strongly. The 
only assumption is that low-energy degrees of freedom 
are fermions that are described by Landau's Fermi liquid 
theory; the extension to, e.g., the superfluid A phase of 
•^He [9] is deferred to a future question, where Nambu- 
Goldstone bosons associated with the spontaneous break- 
ing of U(l) particle number must also be taken into ac- 
count. (In such a Weyl superfluid, the relation between 
the topology of a Fermi surface and anomalies was stud- 
ied theoretically in a different way and verified experi- 
mentally; see Ref. [9] and references therein.) 

As is evident from our arguments, triangle anomalies 
in a Fermi liquid have a "kinematic" origin, independent 
of the details of the Hamiltonian. Namely, we will show 
that Berry curvature modifies the commutation relation 
of the particle number density operator, and that this 
commutator is related to the anomalous Hall effect and 
the triangle anomalies for the fermion numbers near one 
Fermi surface. 

Hamiltonian formulation of Landau 's Fermi liquid the- 
ory. — The fundamental equation of Landau's Fermi liq- 
uid theory is a kinetic equation governing the time evo- 
lution of the occupation number of quasiparticles np(x). 
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where Cp = e^ + Sep, e" is the energy of a single quasipar- 
ticle excitation with energy p, and Sep is the modification 
of its energy due to interactions with other quasiparticles. 
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(3) 



Sriqiy) = ?iq(y) — ?iq IS the deviation from the ground 
state distribution function and /(p,q) are Landau's pa- 
rameter. Above we have neglected the collision term. 

For the purpose of generalizing Landau's Fermi liquid 
theory to systems with Berry curvature, we reformulate 
the kinetic equation as the evolution equation of a Hamil- 
tonian system. In this formulation, the kinetic equation 
has the form 



dtnp{x) = i[H, np(x)], 



(4) 



where the Hamiltonian H is the conserved energy, 

1 /■ dp dq dx 
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and the commutator is postulated as 

d 



/(p, q)(5np(5nq, 



(5) 



d 



[np(x), nq(y)] == -i(27r)3 — <5(p - q) • — ,5(x - y) 

X [^p(y)-"qW]. (6) 

It is straightforward to verify that Eqs. (4), (5), and (6) 
imply Eq. (2). 

The commutation relation (6) is remarkable in the fol- 
lowing respect: assume we have two operators, A and _B, 
linear in occupation numbers, 

A= /"^E^yl(p,x)np(x), ^=/^^B(p,x)np(x), 



(2^)3 
then its commutator will be 



(27r)3 



[AS]=-iy^|^{A, i?}np(x), 
where {A, i?} is the classical Poisson bracket 



(7) 
(8) 



r . . . ^ dA dB dA dB 
R5}(p,x). -.---.-. (9) 

The presence of the Berry curvature, as we shall see, 
changes the classical Poisson bracket and leads to a mod- 
ification of Landau's Fermi liquid theory. 

Berry curvature and Poisson brackets. — Before tack- 
ling the many-body physics of Fermi liquids, let us con- 
sider a single quasiparticle in a theory with Berry cur- 
vature of the Fermi surface. Such a quasiparticle is de- 
scribed by the action [10, 11] 

S = fdt [fx' + A,{x)x' - A,{p)p' - H{p, x)] , (10) 

where H{p,x) is the Hamiltonian whose form is not im- 
portant for us right now, Ai is the electromagnetic vec- 
tor potential, and Ai{p) is a fictitious vector potential 
in momentum space. Combining p and x into one set of 
variables ^°, a ~ 1, . . . , 6, the action can be written as 



S^ dti-LOaiOC-Hm- 



in) 



The equations of motion that follow from this action are 



U^abC = -daH , 



(12) 



where ujab ~ daUJb — dbUJa and da = d/d^"^. We can 
reinterpret this equation as 



r = {^,n = -{r, THi^, 



(13) 



where the Poisson bracket is defined as 



(14) 



where w ^ is the matrix inverse of LOab- For the ac- 
tion (10), the Poisson brackets are [11] 



{Pu Pj} 
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r 1 __ ^ijk^ ^k 



Sij + VtiBj 






(15a) 
(15b) 
(15c) 



where Bi = eijkdAk/dxj, fli = eijkdAk/dpj. 

The invariant phase space is (here ui = detwab) [10] 



dr = V^d^ = (1 + B-rj 



dpdx 



(16) 



It is now clear how to incorporate Berry curvature into 
Landau's Fermi liquid theory. One makes a phase-space 
modification to the Hamiltonian (5), keeps the evolution 
equation (4) unchanged, but alters the commutator of 
np(x) to be consistent with Eqs. (15). We shall now 
work out this commutator. 

Let us assume that there are two operators A and B 
defined as 



d^V^AiOniO, B= d^V^BiOniO. (17) 



A 



Then it seems natural to define the commutator between 
n(^) so that 



[i, B] == -i / d^V^Lj^^daAdbBniO- 



(18) 



This form, however, is deficient in one respect: it makes 
use of n(^) in the whole Fermi volume, while we expect 
the physics to be concentrated near the Fermi surface. 
We shall therefore postulate another form for the com- 
mutator, 

[A, B] = -^^jdi^uj^'^AdaB - BdaA)dbn{i). (19) 

If we integrate by part in this equation, using 
db{\/ujuj°-'^) = (a consequence of the Bianchi identity), 
we bring Eq. (19) into the form of (18). However, now 
the commutator depends only on the physics near the 
Fermi surface. Moreover, we may have problems defin- 
ing the integral in Eq. (18) when the Berry curvature is 
singular inside the Fermi volume (as in the case when the 
Berry curvature flux is nonzero), while Eq. (19) is com- 
pletely well-defined in this case. We will take Eq. (19) as 
the equation defining the commutators of the occupation 
number operator. 

It is possible to write down explicitly the commutator 
[np(x), nq(y)]. We shall not do it here. Instead, we shall 



notice that if A and B are not linear in n(^), but are 
instead general functionals of n(^), then the commutator 
between them can still be computed explicitly, 



|i.Bl.-i/dCv^.-(^8.j^. 



Equation (20) is particularly useful when A is the Haniil- 
tonian, for which we know that 6H/6np{x) ~ ep(x). 

Commutator of density operator. — We now show that 
the Berry curvature leads to an anomalous term in the 
equal-time commutator of the density operator ri(x) at 
two points. Moreover, if the Berry curvature has a 
nonzero magnetic flux through the Fermi sphere, then the 
commutator has a contribution from the external mag- 
netic field. 



[n(x), n(y)] = -i ( V x ^ + — B ) • V(5(x - y), (21) 



where cr is defined as 



"^^"^ = -1(0 



dp „ ^"^p(x) 
(2^)3^^"'= dp, ' 



(22) 



and k is the monopole charge inside the Fermi surface. 



k^— I ds-n. 

2tt 



(23) 



We note that both cr and k involve only the physics near 
the Fermi surface. 

To derive Eq. (21), first we write the density operator 
as 

^(y) =/^ (1 + ^ • ")"p(y) =/dr<5(x - y)np(x). 

(24) 
The commutator of the density operator at two different 
points is, according to Eq. (19), 



/ 1 ^^P 

^^*' ""'^ dx, 



dup 

{X^, Pj}-^ 

dp, 



[n(y), n(z)] = -'-J dV 5{^-y)d^5{^-z) 

-(yf^z). (25) 
The {xi, Xj} term in the commutator is reduced to 
-ia.(5(y-z)y'^e.,fef7fe|^ = -i(Vxa).V<5(y-z), 

where cr is defined in Eq. (22). The {xi, pj} term in the 
commutator can be rewritten as 



dp drip 



■.B.d.5iy-.)j^,n 



' dp J 



(27) 



and, by integration by part, taking into account di^i = 
around the Fermi surface, rip — 1 deep inside the Fermi 
surface and rip = Q far outside the Fermi sphere, it be- 
comes 



-iAB.V<5(y- 



(28) 



Combining two contributions, we find Eq. (21). 

From density- density commutator to anomalous non- 
conservation of current. — The connection between the 
anomalous density-density commutator [the term propor- 
tional to B in Eq. (21)] and triangle anomalies is known 
in the context of relativistic quantum field theory [12, 13]. 
Here we derive this connection using the Hamiltonian for- 
malism and show how the anomalous Hall current and the 
triangle anomaly can be traced to the two contributions 
to the density-density commutator. 

Let us first assume that our system is in a static mag- 
netic field, but the electric field is turned off. In this 
case, the system is described by the Hamiltonian (5), and 
by commuting the Hamiltonian with the particle number 
operator n(x), the continuity equation can be derived, 



i[iJ, n] 



Vj, 



(29) 



where the particle number current j is 



dp 

(2^)= 



dup 
'P9p 
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ax 



Note that by integration by part, the first term in the 
brackets in the right-hand side of Eq. (30) can be written 
in the familiar form ripV, where v = dep/dp. This would 
be the only term in the current in the absence of Berry 
curvature. 

Now we turn on a static electric field by putting the 
system in an external scalar potential 0(x), E = — Vi/). 
The Hamiltonian is now 



H' = H 



dx(/)(x)n(x). 



(31) 



The added term docs not commute with n and changes 
the time evolution of the latter. 



dtn{x) ^ i[H', n(x)] = -V-j- iVxa + —^B \ • V(/)(x). 

(32) 
This equation can be rewritten as 

(33) 



dtn + V-j' = ^E-B, 



where 



j'=j + Excr. 



(34) 



The second term in Eq. (34) is the usual anomalous Hall 
current. On the other hand, Eq. (33) implies that the 



particle number around the Fermi surface is not con- 
served when both electric and magnetic fields are turned 
on. This is the effect of triangle anomalies in quantum 
field theory. For example, relativistic right-handed free 
fcrmions have k — 1, and left-handed free fermions have 
fc = — 1. Here we show that this effect depends only on 
the monopole charge of the Berry curvature on the Fermi 
surface and is not modified by interactions. Since the to- 
tal charge is conserved, all different contributions to the 
current nonconservation should sum up to zero. 

Chiral magnetic effect — Let us compute the current, 
given by Eq. (30) , in the thermal equilibrium state, where 
quasiparticles have a Fermi-Dirac distribution function. 



/(^) 



1 



£p-M 
T 



(35) 



There are three contributions to j corresponding to 
three terms in the brackets in the right-hand side of 
Eq. (30). The third term involves spatial derivatives and 
vanishes in the ground state. We now show that the first 
term also vanishes identically. For this end it is useful 
to introduce the function g(x) — J_ dyyf'(y), for which 
g{—(yo) = g(+oo) = 0. Then 



dp 



drir, 



(2^)3 ^P dp 



dp 

{2nY 
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Similarly, the second contribution can be written as 



B 



dp 



n 



"If-^l^'^) 



0. 
(36) 

(37) 



and the integrals can be evaluated as in Eq. (27). As the 
result, we find 



J - 4^^^- 



(38) 



Let us assume for dcfinitcncss that there are two Fermi 
surfaces with fc = 1 and fc = — 1. If the two Fermi sur- 
faces have equal chemical potential, then the total cur- 
rent is equal to 0. However, if the chemical potentials are 
unequal (which can be achieved by turning on an E • B 
for a finite time), then there will be a current equal to 
(;ti+ — /L(_)B/47r^ in this quasiequilibrium state. This is 
the chiral magnetic effect [14, 15]. 

Conclusion. — The calculation above ties the anoma- 
lous current nonconservation to a property of the Fermi 
surface only (the Berry curvature) and hence can be ap- 
plied to Fermi liquids, even when the interaction between 
original fermions is strong. This is done by using a kinetic 
equation for quasiparticles; a more microscopic deriva- 
tion of this equation [16, 17] is desirable. 

It would be interesting to explore further physical con- 
sequences of Berry curvature on Landau Fermi liquid the- 
ory. Particularly interesting are the effects of Berry cur- 
vature on the collective modes and on the response of the 



Fermi liquids. It is also interesting to include the collision 
term into the kinetic equation and investigate the hy- 
drodynamic regime. In relativistically invariant theories, 
the effects of triangle anomalies have been investigated, 
both within hydrodynamics and by using gauge-gravity 
duality [18-20], and there have been attempts to derive 
hydrodynamics from kinetic theory [21]. The kinetic ap- 
proach allows us to go beyond the hydrodynamic regime 
and beyond systems with relativistic invariance. On the 
other hand, some interesting phenomena associated with 
anomalies in relativistic theories, like the Alfven-type 
modes propagating along the direction of the magnetic 
field [22, 23], may be directly investigated using the ki- 
netic equation in the condensed- matter context. 

Finally, the understanding obtained here should allow 
one to formulate the criteria of anomalies matching for 
dense states matter, i.e., quark matter phases with Fermi 
surfaces. 
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